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Lattice results are available for AS = 2 matrix elements for the first time in full QCD, which 
improve considerably the status of hadronic uncertainties in K-K mixing with respect to earlier 
phenomenological studies. Using an average of the ETMC and RBC-UKQCD results, we analyze 
ck in Natural SUSY. This scenario arises as a consistent BSM framework after the latest results 
from the LHC. The analysis is improved with respect to previous studies including next-to- leading 
order matching conditions of order o? s . We derive new bounds for SUSY mass insertions in the 
scenario with a light third generation and study the implications for squark and gluino masses, 
compared with direct searches at the LHC. Assuming natural values for the flavor violating SUSY 
couplings of both chiralities, we find that the sbottom must be heavier than 3 TeV for a gluino mass 
up to 10 TeV. In this scenario no natural values for squark and gluino masses can satisfy the flavor 
bounds. 



1. INTRODUCTION 

Flavor physics observables related to mixing and de- 
cay of K, D and B mesons pose strong bounds on New 
Physics (NP) models. A strong constraint on the scale of 
NP comes from the measurement of 6k, related to indi- 
rect CP violation in the neutral Kaon system, which sets 
a lower bound on the NP scale around A ~ 10 4 TeV in 
the presence of flavor violating couplings of 0(1) [3, Q. 

In order to study flavour observables, one has to face 
the calculation of the matrix elements of the relevant 
local operators. In the case of ck, the matrix elements of 
the full set of AS = 2 operators beyond the SM has been 
recently computed in full QCD by the ETMC and RBC- 
UKQCD lattice QCD collaborations 0,3. These results 
constitute a considerable improvement with respect to 
previous results in the quenched approximation [5| . The 
model-independent bounds on the scale of New Physics 
imposed predominantly by the operator Q4 has increased 
almost by a factor of 3 [3[. 

These new results can be immediately used to set con- 
straints on Supersymmctry, putting bounds cither on its 
flavor violating couplings, or on the SUSY masses, if some 
particular scenario is chosen for the flavor violation. A 
first analysis has been performed in Ref. @, where they 
consider the QCD running between the scale set by the 
heavy squark masses, and the lower scale set by the gluino 
mass (and eventually a light third generation). 

After the first run of the LHC, direct SUSY searches 
have established relatively strong bounds on the masses 
of the squarks of the first two generations, more moder- 
ate bounds on the gluino mass, and still weaker bounds 
on third generation squarks. This circumstance is in fair 
connection to the spirit of Natural SUSY, where the only 
strongly interacting SUSY partners required to be light 
are the squarks of the third generation, and to a lesser 
extent, the gluino. This generic SUSY scenario is con- 
sistent with naturalness and with current results from 
direct searches at the LHC 0-0- 



In this letter, we study the bounds imposed by ex on 
Natural SUSY taking into account the recent lattice QCD 
results for the matrix elements, as well as NLO match- 
ing conditions for the AS = 2 Wilson coefficients. We 
begin in Section [2] reviewing briefly the relevant formulae 
for ex beyond the SM, and in Section [3] we combine the 
two different sets of lattice QCD results for the matrix 
elements, obtaining averaged results to be used in the 
phenomenological analysis. In Section [4] we summarize 
the relevant details concerning flavor violation in Natu- 
ral SUSY, in Section [3] we study the constraints on the 
flavor violating couplings and in Section [5] we study the 
implications on squark and gluino masses, under certain 
generic assumptions concerning the flavor violation. 



2. KAON MIXING IN THE SM AND BEYOND 



The parameter ex is given by fioj 
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where is the short distance contribution at the 

charm scale. Assuming non-relativistic normalization for 



matrix elements, M$ = (K°\H eS \K°). This short dis- 
tance contribution can be split into the SM and NP com- 
ponents: 
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where the NP contribution can be related to the SM and 
experimental values for e^- through 



ImM NP = y^Arr^ 



(|eA-r P -|eA-| bM ) 



(3) 



with l/ Ke = l/(V2sin0 e )(l - (>£/(V2\e% p \)), namely, 
k £ = 0.94 ± 0.02 H3]. For the SM value of e K we take 



SM 



(1.9 ±0.3) x 10" 



(4) 
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computed in Ref. [TTJ and rescaled to our value of k c and 
to the more recent average of the B parameter Bk given 
in Ref. [3 (see Section©. 

The current experimental values |eif| exp and Am^ p 
are given bv[l3j 

Am^ p = (3.483 ± 0.006) • 10~ 15 GeV , (5) 
\e K \ cxp = (2.229 ± 0.010) • 10" 3 , (6) 

which together with Eq. @ imply the following bound 
on the NP contribution: 



ImMf/ = (1.7 ± 1.6) • 10~ iS GeV 



(7) 



The most general effective Hamiltonian for K-K mix- 
ing beyond the SM is given by 



n cB = J2 C * (aO 2* (m) + E c * M Qi (m) . ( 8 ) 

i=l i=l 

where the SUSY basis of operators is 

Qi = d a ~f^P L s a dp^P L sp 

Q 2 = d a PLS a dpP L Sfj 

Q3 = d a P L s^dfiP L s a (9) 

Qi = d a P L s a dpP R sp 

Qb = d a P L spdpP R s a 

together with the chirally-flippcd operators Qi,2,3 ob- 
tained from Qi.2,3 with the substitution L <-> R. The 
chiral projectors are defined as Pl,r = (1 T 7s)/2. 
The New Physics amplitude is then given by 



< 2 p = E^ p (^)(^i^(m)i^ ) 



(10) 



The matrix element for the SM operator Qi is related 
to the bag parameter Bk (in the non-relativistic conven- 
tion) 



(K°\QiOj,)\K°) = \m K f K B K {p) , (11) 

and the matrix elements of the operators Q2,3,4,5 are usu- 
ally normalized to (Qi), defining the ratios Ri as 



e?\Qi(n)\K°) 

"iM = -=0; 



(12) 



The NP Wilson coefficients must be given in the same 
renormalization scheme as the matrix elements, and at 
the same renormalization scale [i. Since the matching 
conditions are computed at the matching scale A related 
to the masses of the heavy particles, the Wilson coeffi- 
cients must be evolved down by means of the Renormal- 
ization Group. The evolution matrix at NLO in QCD for 



A = 1 TeV - RI scheme 


61(A) =0.762 


62(A) = 2.544 


6s(A) = -0.002 


62(A) = -0.591 


63(A) =0.390 


64(A) = 4.823 


6s(A) =0.186 


64(A) = 1.351 


6s(A) = 0.875 



TABLE I: Values from ref. [3 for the NLO AF = 2 evolution 
coefficients from A = 1 TeV to [i = 2 GeV in the Landau RI 
scheme and in the SUSY basis. 



fj, = 2 GeV in the RI scheme is given in Ref. [14[ . Taking 
this into account, we can write 



M NP 

iu 12 



\m K f K B K [ai(A)(Cf p (A) + C7HA)) 

+ [6a(A)ifc + 6 3 (A)i? 3 ](C 2 NP (A) + C 2 NP (A)) 
+ [6a(A)ifc + 6 3 (A)i? 3 ](C 3 NP (A) + C 3 NP (A)) 
+ [64(A)i?4 + 6- 5 (A)i? 5 ]C 4 NP (A) 



[64(A)i?4 + e 55 (A)i? 5 ]C 5 NP (A) 



(13) 



where the NLO evolution coefficients £y (A) for A = 1 
TeV in the Landau- RI scheme are collected in Table Q] 
Including NLO matching conditions for the Wilson coef- 
ficients C{, the combination 6j(A)C7j(A) is independent 
of the matching scale at NLO. Eqs. Q and (jTB")) will be 
used in the following sections to study the constraints 
from €k on NP. 



3. REVIEW OF LATTICE QCD RESULTS FOR 
AS = 2 MATRIX ELEMENTS 

The bag parameter Bk has been calculated in full 
QCD by lattice groups since 2004 [l5[. The average re- 
sult up to 2010 for the corresponding renormalization- 
independent parameter Bk is given by fl6| : Bk = 
0.738(20). Recently, new refined lattice studies have be- 
come available I U7H20| . Here, we use the updated world 
average of Ref. [12| : 



B K = 0.7643(97) 



(14) 



This leads to the following value for the B-parameter in 
the Landau-RI renormalization scheme: 



B^\2 GeV) 



0.546(7) 



(15) 



This year, the ratios Ri in Eq. (ri2|) have been calcu- 
lated in full QCD for the first time, by the ETMC and 
RBC-UKQCD collaborations H H, with N f = 2 and 
2+1 active flavors respectively. These results supersede 
previous ones in the quenched approximation |2ll . [22| . 

The RBC-UKQCD and ETMC matrix elements are 
given in the SUSY basis, at a renormalization scale of 
3 GeV and in the MS scheme of Ref. [24j ]. We perform 



3 



MS at fi = 3 GeV 




ETMC [3] 


RBC-UKQCD [4J 


Our Average 




-16.3(0.6) 


-15.3(1.7) 


-16.2(0.9) 


Ml*) 


5.5(0.4) 


5.4(0.6) 


5.5(0.3) 


Ra(h) 


30.6(1.3) 


29.3(2.9) 


30.4(1.2) 


Rsif*) 


8.2(0.5) 


6.6(0.9) 


7.8(1.2) 



the chiral basis and the MS scheme of Rcf. 24 



TABLE II: Unquenched lattice QCD results for the ratios 
Ri(fi) as given in Refs. @, Q|, in the MS scheme of Ref. 
given at the renormalization scale fi = 3 GeV. Our average 
results are computed as explained in the text. 



a weighted average of both results using the procedure 
described in Ref. [l3[. In the case of R2 and R$ the 
RBC-UKQCD central values lie outside the average error 
band. To take this into account we increase the errors 
to include the RBC-UKQCD central values. The RBC- 
UKQCD and ETMC results together with our averaged 
values are collected in Table [TTI and represented in Fig.Q] 

For the phenomcnological analysis — according to 
Eq. (fTB"]) — one has to combine the Ri averages in Table [TTI 
with the £ij(A) factors in Table HJ For this purpose, how- 
ever, the quantities Ri and £ij(A) must be defined using 
the same renormalization prescription. The ratios Ri in 
Table |TTJ are defined in the MS scheme of Ref. [U] and 
at [i = 3 GeV, while the coefficients ^ (A) in Table UJ 
are given in the Landau-RI scheme for \i = 2 GeV. We 
find it more convenient to transform the ratios Ri to the 
Landau-RI scheme at /1 = 2 GeV, e.g. to the prescription 
in which the coefficients £y (A) are given. 

We first perform the QCD running from 3 GeV down 
to 2 GeV. This running is performed at NLO by means of 
the two-loop anomalous dimensions given in Ref. (24| (see 
also Ref. [2q|). However, in Rcf. [24| the renormalization 
is carried out in the so-called chiral basis of operators, Qi 
(see Eq. (2.1) of Rcf. [24|). The translation between both 
bases is a fierz transformation: Qf ch (/i) = A»jQ*- (/i), 
where the matrix A is given by 

/ 1 \ 

0-2 

A= 1 . (16) 

10 

\0 4 8 / 

where the label 'sch' stands for either scheme, MS or RI. 
Fierz transformations introduce a different prescription 
for evanescent operators in the MS scheme, which makes 
the MS scheme of Ref. [j| used by RBC-UKQCD and 
ETMC different from the MS scheme in Refs. [25|,|40[. 

The QCD running from 3 GeV down to 2 GeV is 
given by Qf* (2 GeV) = 17(3 GeV, 2 GeV) ji Qp(3 GeV), 
where U(3 GeV, 2 GeV) is the NLO evolution matrix in 



17(3,2) 



f 1.035 

1.022 0.011 

0.130 0.830 

0.887 -0.474 

V 0.001 1.152 / 



, given by 
\ 

• (17) 



The evolution is performed in 4-flavor QCD [26[, consis- 
tent with the fact that the charm quark is a dynamical 
degree of freedom from 3 to 2 GeV for the NP contribu- 
tions parametrized in Eq. (| L3[) . The value of the strong 
coupling at these scales is obtained from a s (m c ) running 
up to 2 and 3 GeV in the 4-flavor theory. We use the 
full results for the running of a s (/x) from Ref. [23[ , giving 



0.3041 and ai 4) (3GeV) = 0.2552. The 



ai 4) (2GcV) 

relevant inputs at the charm scale are a s (m c ) = 0.3537 
and m c (m c ) = 1.28(1) GeV [27j]. 

The conversion to the RI scheme is performed by 
means of the NLO matrix Ar^s_^ RI of Ref. [23 |. namely 

Qf r (/i) = Ma QpG*), with M = [1 — (a s /47r)ArM^ RI ] 
and ArHs^Ri given by 



Arv 



/ 0.879 




V 





-1.129 
0.307 









-6.773 
10.871 







\ 



(1 

5.644 0.214 

12.939 2.689 / 



This matrix can be rotated to the SUSY basis by means 
of the rotation A given in Eq. (fTBf . The result will differ 
from the one in Refs. [25|,[4(j because the MS renormal- 
ization scheme is not the same. 

Summarizing, to work out the ratios Ri(fi) at /i = 2 
GeV in the Landau-RI scheme from Ri(fi) at fi = 3 GeV 
in MS we make 



R™ (2 GeV) = i?j MS) (3 GeV) 
where the transformation matrix Mu is defined as 



Mi 



[A- 1 MU T (3,2)A], 
[A- 1 MC/ T (3,2)A] ] 



(19) 



(20) 



Numerically, we find: 



/l 



M 







Vo 















\ 



0.743 -0.037 

0.073 1.083 

0.608 -0.001 

-0.131 1.037 ) 



(21) 



Applying this transformation to the averaged lattice re- 
sults of Table [TTI we get 



4 RI) (2 GeV) 



p(R-i) 
i?< RI> 



(2 GeV) 
(2 GeV) 
(2 GeV) 



-12.2(0.7) , 
4.8(0.3) , 
18.5(0.7) , 
4.1(1.2) . 



(22) 



4 



RBC-UKQCD'12 
R 2 = -15.3(17) 



ETMC'12 
R 2 = -16.3(06) 



-16 



-14 



-12 



-10 



RBC-UKQCD'12 
R 3 = 5.4(06) 



ETMC'12 
R 3 = 5.5(04) 



-2 



4 6 



10 12 



RBC-UKQCD'12 
tf 4 = 29.3(29) 



ETMC'12 
R 4 = 30.6(13) 



RBC-UKQCD'12 
R 5 = 6.6(09) 



ETMC'12 
R 5 = 8.2(05) 



26 28 30 32 34 36 38 2 4 6 8 10 

FIG. 1: Weighted average of unquenched lattice results of Refs [3, Q, for the ratios of matrix elements Ri, given in the MS 
scheme of Ref. at the renormalization scale fj, = 3 GeV. The lighter colored error bars correspond to the enlarged errors 
that include the central values of the RBC-UKQCD results. See the text for details. 



These values, together with B { P(2GeV) = 0.546(7) of 
Eq. (|15[) . will be used in the phenomenological analysis 
in Sections [S] and El 

4. FLAVOR VIOLATION IN NATURAL SUSY 

In general SUSY models, flavor violation in the quark 
sector is mediated predominantly by strong interactions, 
via flavor-changing quark-squark-gluino interactions in- 
duced by soft SUSY-breaking terms, flavor physics. 

Let M. q be the squark mass matrix in the q = u,d 
sector, given in the super-CKM basis. In order to go 
to a physical basis where squarks do not mix with each 
other, a rotation is performed in the squark sector alone 
to diagonalize the squark mass matrix: 

m 2 diag =r q M 2 q rl, (23) 

After this rotation is performed, the 6x6 unitary matrix 
T q appears in the quark-squark- gluino vertex: 

£ g g S = -\/2<? s rf (q J T a q l )g + h.c. (24) 

where q l are three left-handed (i=l,2,3) and three right- 
handed quarks (i=4,5,6), of type q = u or d. This vertex 
leads to squark-gluino loop penguin and box diagrams 



that contribute (among other things) to AF = 1 and 
AF = 2 processes. As an example, the contribution 
to a s — > d transition is given at the leading order by 
A™Z ~ (a./m s )I*/*r£/(ro?/ro§), where f(x) is a pen- 
guin function. It is clear that both in the case of degen- 
eracy (fhi = rh) and in the case of alignment (r 1 - 7 = <5.y ) , 
the amplitude vanishes. 

Mechanisms suppressing flavor violation in SUSY such 
as degeneracy (SUSY-GIM mechanism) or alignment are 
required by flavor physics data, if the soft SUSY-breaking 
scale is low to comply with naturalness. In the absence 
of such mechanisms, the NP scale must be as high as 
A > 10 4 TeV in order to satisfy bounds from K — K 
mixing P, Q (the bounds from _B-physics are some- 
what weaker A > 10 2 TeV). The absence of a natural 
symmetry-based principle providing a sufficiently effec- 
tive suppression of flavor violation in the presence of a 
low SUSY scale, without challenging naturalness, is a 
manifestation of the SUSY flavor problem. 

However, naturalness does not require all the soft 
masses to be low, but only those linked more strongly 
to the Higgs. In the strong sector, the stops Il,r con- 
tribute at one loop to the higgs mass and should be not 
much heavier than about ~ 500 GeV, while the gluino 
contributes at the two-loop level and should not be heav- 
ier than about ~ 1.5 TeV @, H[i assuming that the fine 



5 



tunning is not worse than ~ 10%. By SU{2)l symme- 
try, the "left-handed" sbottom b\ is also required to be 
light. Beyond these restrictions, first and second gener- 
ation squarks can be heavy, providing a scale suppres- 
sion to flavor violation without compromising natural- 
ness. These type of SUSY models have been collected 
under the name of Natural SUSY. 

In Natural SUSY, the transition sod mediated by 
first and second generation squarks is suppressed by their 
heavy masses, and the competing process where the tran- 
sition is mediated by third generation squarks takes over, 
even though it is second order in flavor violating cou- 
plings. This mechanism relates flavor violation in K and 
B physics. K — K mixing sets bounds on flavor violating 
couplings related to the third family, that are compara- 
ble to those derived from B physics [28|, H3- However, 
for this mechanism to work, the scale suppression pro- 
vided by the squark masses of the first two generations 
is in general not enough, and an additional U (2) flavor 
symmetry might be invoked [29j . 

Taking into account these considerations, we consider 
a Natural SUSY scenario with first generation squarks 
of mass ~ fhh around ~ 10 TeV, and third generation 
squarks of mass ~ fh( around ~ 500 GeV. A suitable pa- 



rameterization of the rotation matrices T q = 



is given by [3(J, |31| 
/ 



1 
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5 LL °e 



r(/,23* 

s ee 



_<5?,13 
LL 

_?9,23 
LL 

C<9 






-s e e 



-itp 



°rr s ee 
1 


g,13* 



°rr b ee 

1 

q,23* 






see 



.$9,13 
RR 

°RR 



(25) 



eg, 1,5* c-Q-Zj* 

V 6 rr c e °rr c e c e ) 

where cq = cos 9 q and so = sin 8 q , with 8 q the mixing 
angle in the q LR sector. The mass insertions &Yl rr are 
the couplings responsible for the flavor transitions, and 
can be bounded imposing flavor constraints. A similar 
parameterization for rotation matrices with non- degen- 
erate squarks has been considered, for example, in the 
phenomenological analyses in Refs. [32Tl34j . an impor- 
tant difference being that 8f^ L , 8 RR were set to zero to 
kill effects in kaon physics. 

In the next section we consider the bounds that can 
be derived from ex assuming a squark spectrum of the 
type discussed above. On the other hand, these mass 
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FIG. 2: \e K \ in SM + SUSY vs the matching scale for a set of 
SUSY parameters consistent with experimental bounds. The 
blue dashed line is LO and the red solid line is NLO. The 
point at which the LO and the NLO coincide depends on the 
point in SUSY parameter space. 



insertions receive contributions from soft SUSY-brcaking 
parameters in the Lagrangian, as well as from Yukawa 
couplings. Assuming no particular cancellation between 
these two (in principle unrelated) contributions, leads to 
a natural size of the mass insertions that can be used to 
infer bounds on squark and gluino masses. This is the 
target of Section [6] 

In order to study the constraints from flavor observ- 
ables, the SUSY amplitudes must be computed. The 
model- dependent part of these amplitudes is encoded in 
the matching conditions, that is, the values of the Wilson 
coefficients in the effective Hamiltonian at the matching 
scale A. These matching conditions are known to NLO 
in strong interactions: leading order matching conditions 
can be found in Refs. [35| and [36| for AF = 1 and 
AF = 2 processes respectively. Two-loop NLO correc- 
tions to AF = 1 have been computed in Refs. (37l. f38j] . 
while the full NLO corrections to AF = 2 can be found 
in Refs. [Hill. 

While it can be argued that NLO corrections are nu- 
merically small and have no real impact on the bounds 
derived for the SUSY parameters, it should be noted that 
at leading order the amplitude suffers from a substantial 
renormalization scale dependence that leads to large un- 
certainties. The two main reasons for this sensitivity to 
the renormalization scale are [1^, H3] : (a) the leading or- 
der contribution is proportional to a 2 s , while there is no 
definition of the renormalization point at LO, and (b) 
the anomalous dimensions of the operators in Eq. (fT0|) 
are large. 

In order to stress this point we show, in Fig. [21 the de- 
pendence of \ck \ on the SUSY matching scale comparing 
the LO and NLO results. There is clearly a considerable 
reduction in the renormalization scale ambiguity when 
going from a LO to a NLO matching. By performing a 
complete NLO analysis, it is justified to ignore the un- 
certainty related to the variation of the renormalization 
scale. We emphasize that a complete NLO analysis in 
non-degenerate SUSY scenarios has never been done be- 
fore, and we also note that, in general, existing LO anal- 
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yses do not take into account the renormalization scale 
uncertainty. 

Besides the renormalization of a s and squark and 
gluino masses that must be taken into account at NLO, 
flavor changing renormalization of quark and squark 
propagators have to be considered. The (finite) renor- 
malization of quark fields induced by squark-gluino loops 
leads to chirally-enhanced effects that can be numerically 
important (see Refs. [ij],[42[). However, in an "on-shell" 
scheme for the super-CKM basis these corrections are 
absent. The difference between both schemes boils down 
to a different definition for the mass insertions (see ap- 
pendix C of Ref . [39[ ) . In this letter all mass insertions 
are defined in the on- shell scheme. The (infinite) renor- 
malization of squark fields induced by the squark tadpole 
implies that the diagonalization of the squark mass ma- 
trices must be performed at each renormalization scale. 
We therefore define the rotation matrices T q (ft,) at a fixed 
scale ft, and include in the matching conditions the con- 
tribution from non-diagonal squark masses, which are of 
order my(/i) ~ a s \og ft/ (i. These in fact contribute to 
the RG equation and to the reduction of the renormal- 
ization scale uncertainty 

Apart from strong-interaction squark-gluino correc- 
tions, contributions from chargino-squark loops are rele- 
vant in certain scenarios due to the role of A-terms. We 
understand that both contributions are mostly uncorre- 
lated in a general set-up, meaning that both contribu- 
tions set independent bounds on SUSY (see for example 



0.10 



Section 3 of Ref. |43[). In this letter we focus on the 



conclusions that can be taken from squark-gluino contri- 
butions alone. A study of the effect of chargino contri- 
butions is certainly worthwhile, but beyond the scope of 
this note. 



5. CONSTRAINTS FROM e K ON FLAVOR 
VIOLATING COUPLINGS 

In this section we derive constraints on the insertions 
^ll rr an d $ll rr from the measurement of zk- The 
bounds are obtained imposing the constraint in Eq. Q 
on the NP amplitude of Eq. (fT3")) . where the NLO match- 
ing conditions for the coefficients C;(A) are taken from 
Ref. [39| . The matching scale is fixed at A = 1 TeV, 
which is justified at NLO according to the discussion 
in the previous section. The coefficients C, depend on 
the gluino mass m§, the heavy and light squark masses 
to/j, fhi, and the rotation matrices T q , all defined at the 
matching scale. For the rotation matrices we use the 
parameterization of Eq. (j2"5)l . 

For the analysis we fix the masses to m/, = 10 TeV, 
m bL = 500 GeV, fh bR = 700 GcV and m s = 1 TeV. We 



put the flavor violation in the up sector to zero, and we 
consider two scenarios: 5 RR — (LL only) and 8^ L = 
6t R (LL=RR). 

In Fig. [3] (upper plot) we show the one and two-sigma 
constraints on the 8f/ L — 5f/ L plane in the case of LL 
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FIG. 3: Constraints on mass insertions from \ex\ in the case of 
LL mixing only (upper plot) and LL=RR (lower plot), at one 
sigma (dark red) and two sigma (light red). The constraints 
are obtained as explained in the text. 



mixing only. The constraints are obtained by a standard 
X-square minimization, in the situation where the com- 
plex phases make the amplitude purely imaginary. These 
constraints are therefore the most conservative on the 
magnitude of the mass insertions, \6^ L \. These bounds 
can be approximately summarized by the constraint: 



Im[(5 



db jrsb** 
LL°LL, 



2 ] < 1.7 - 10~ 6 at 95%C.L. 



(26) 



The one and two sigma constraints in the case of LL=RR 
mixing are shown in the lower plot of Fig. [3] In this 
case, the relevant bound can be expressed approximately 
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Im [6?J%6$ R 5j&] < 1.6 • 10- 9 at 95% C.L. (27) 

These approximate results are obtained neglecting terms 
in the amplitude containing a product of more than four 
mass insertions. Since the mass insertions are small, and 
having checked that the numerical coefficients of such 
terms are also small, this approximation is fully justified. 



6. IMPLICATIONS FOR SQUARK AND GLUINO 
MASSES 

Focusing on the LL sector, the squark mass matrices in 
the super-CKM basis are given by M\ LL = VurfigV^ and 

M 2 d LL = VdfriQVj. Here, are the soft masses for the 
squark SU{2)l doublets, and the matrices V u .d are the 
rotations transforming left-handed quark supermultiplets 
from the weak to the super-CKM basis, and such that 



VuVl 



Vckm is the CKM matrix. The link between 
M^ L and M^ L imposed by SU(2)l symmetry is then 
M\ LL = Vqkm M\ LL Vj KM . This relationship has been 
used to relate flavor violation in K— K and D — D mixing 
(see for example Refs. 0, 5EJ ) . 
We can diagonalize both matrices applying the rota- 



tion matrices T 



(LL) 



in the LL sector: 



m 2 



diag 



t ~ 2 

r„ v C km r d m diag 



1 d CKM 1 u 



(28) 



of order 



'diag 



= diag^m^m^ifnf) up to perhaps terms 
We note that we have dropped the su- 



perscript {LL) in the T q matrices. For convenience, we 
define W = T u V C km T\. The right hand side of Eq. ([28]) 
can be written as 



(29) 



r.h.s 



Expanding in the same way the left hand side, Eq. ([28]) 
leads to Uzi — 5^ + 0(rh^/m^). This equation sets a 
natural size for the mass insertions. For example, in 
the case in which the up quark and squark sectors are 
approximately aligned, we have T u ~ 1 and therefore 
0(m 2 /fh\), which translates into [3(| 



r 3 

1 d 



"CKM 



$LL ~ V CKM + ^(^MD 



(30) 



This discussion is justified when the ratio mj/rn^ is 
very small. On more general grounds, the condition that 
any chiral-conserving entry of the matrix M. q is at least 

of size rhf, leads to III > rhf/fh 2 , H3l- Excluding unex- 
pected cancellations, we expect 



°LL ~ max ( V CKM> m ^/ m feJ 



(31) 



In this section we assume that the mass insertions sat- 
isfy the lower bounds S L '^ RR > V^icM' an< ^ study the 



> 

o 
H 



I 0.3 



> 
H 




0.0 L_ 

0.0 0.1 0.2 0.3 0.4 0.5 Of 

Mgiumo (TeV) 




M 



gluino 



(TeV) 



FIG. 4: Constraints on the masses of the gluino and the sbot- 
tom (rht = rhb L = rhb R ) from \ek\ in the case of LL mixing 
only (up) and LL=RR (down), at 65% (dark red) and 95% 
C.L. (light red). The dashed line correspond to the 95% C.L. 
constraint obtained from LO matching conditions. The con- 
straints are obtained as explained in the text. 



implications of the measurement of ck on squark and 
gluino masses. 

The results are shown in Fig.|3]in the case of LL mixing 
only (upper panel) and LL=RR mixing (lower panel), for 
heavy squarks of 10 TeV and rhb L = rhb R . Also shown 
are the LO constraints, that turn out to be less stringent 
than the NLO ones. In the absence of RR mixing, for 
a gluino heavier than 200 GeV, the sbottom mass is un- 
constrained. These bounds do not compete with direct 
searches at the LHC. The situation is quite different in 
the case of LL=RR mixing (with rhb L = mj, s ). In this 
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case the operator Q4 gives a big contribution to ex be- 
cause of the chiral enhancement of its matrix element, 
its large anomalous dimension, and because the coeffi- 
cient C4 is numerically large. We find that the sbottom 
must be generically heavier than about 3 TeV indepen- 
dently of the gluino mass (for m- g < 10 TeV). This situa- 
tion is clearly excluded by naturalness. This is an exam- 
ple where the flavor bounds are far more stringent than 
the direct searches at the LHC. An intermediate scenario 
with < 8rr < Sll will lead to constraints that lay in 
between the two extreme situations considered. 



Concerning the SUSY analysis, we show the impact 
of the inclusion of NLO matching conditions, reducing 
considerably the renormalization scale uncertainty. The 
bounds on the flavor violating couplings are summarized 
in Fig. [3] They can be approximated by the bounds given 
in Eqs. (|26[) and (|27[) . Assuming a natural size for mass 
insertions (see Eq. l30|) . we derive lower bounds on squark 
and gluino masses. In the case of LL and RR mixing, the 
bounds are much stronger than the direct bounds from 
the LHC, implying a sbottom heavier than 3 TeV in this 
scenario. 



CONCLUSIONS 

We have analyzed the impact of the latest lattice QCD 
results for AS = 2 matrix elements in full QCD on Natu- 
ral SUSY, with NLO matching conditions for the Wilson 
coefficients. The weighted average of the ETMC and 
RBC-UKQCD results for the matrix elements at 2 GeV 
in the Landau- RI scheme are collected in Eq. (|23|) . They 
imply a big progress compared to older quenched results, 
and can be used to set constraints on New Physics. 
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